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We present an improved prescription for matter power spectrum in redshift space taking a proper 
account of both the non-hnear gravitational clustering and redshift distortion, which are of particular 
importance for accurately modeling baryon acoustic oscillations (BAOs). Contrary to the models 
of redshift distortion phenomenologically introduced but frequently used in the literature, the new 
model includes the corrections arising from the non-linear coupling between the density and velocity 
fields associated with two competitive effects of redshift distortion, i.e.. Kaiser and Finger-of-God 
effects. Based on the improved treatment of perturbation theory for gravitational clustering, we 
compare our model predictions with monopole and quadrupole power spectra of N-body simulations, 
and an excellent agreement is achieved over the scales of BAOs. Potential impacts on constraining 
dark energy and modified gravity from the redshift-space power spectrum are also investigated 
based on the Fisher-matrix formalism. We find that the existing phenomenological models of redshift 
distortion produce a systematic error on measurements of the angular diameter distance and Hubble 
parameter by 1 ~ 2%, and the growth rate parameter by ~ 5%, which would become non-negligible 
for future galaxy surveys. Correctly modeling redshift distortion is thus essential, and the new 
prescription of redshift-space power spectrum including the non-linear corrections can be used as 
an accurate theoretical template for anisotropic BAOs. 



I. INTRODUCTION 

Galaxy redshift surveys via the spectroscopic measure- 
ments of individual galaxies provide a three-dimensional 
map of galaxy distribution, which includes valuable cos- 
mological information on structure formation of the Uni- 
verse. The observed galaxy distribution is, however, ap- 
parently distorted due to the peculiar velocity of galaxies 
that systematically affects the redshift determination of 
each galaxy. The anisotropy caused by peculiar veloci- 
ties is referred to as the redshift distortion, which com- 
plicates the interpretation of the galaxy clustering data 
(e.g., Bi). 

Nevertheless, redshift distortion provides a unique way 
to measure the growth rate of structure formation, which 
has been previously used for determining the density pa- 
rameters of the Universe (e.g., 0, 3), and is now recog- 
nized with great interest as a powerful tool for testing 
gravity on cosmological scales (e.g., [5|-|1|). Redshift dis- 
tortion also provides a helpful information on the dark- 
sector interactions [loj , where the dark e nerg y is dynam- 
ically coupled with dark matter (e.g., (TTI.Il2|). Note that 
the distortion of the galaxy clustering pattern also arises 
from the apparent mismatch of the underlying cosmol- 
ogy when we convert the redshift and angular position of 
each galaxy to the comoving radial and transverse dis- 
tances. This is known as Alcock-Paczynski effect [l3| . 
and with the baryon acoustic oscillations (BAOs) as a 
robust standard ruler, it can be utilized for a simulta- 



neous measurement of the Hubble parameter H{z) and 
angular diameter distance Da{z) of distant galaxies at 
redshift z (e.g., iM-il8|]). 

In these respects, anisotropic clustering data from 
galaxy redshift surveys serve as a dual cosmological probe 
of the cosmic expansion and the gravity on cosmological 
scales, from which we can address properties of both the 
dark energy and modification of gravity responsible for 
the late-time cosmic acceleration. Although current data 
are not yet sensitive enough to sep arately measure H{z), 
Da{z) and growth rate (see 0,millj for current status), 
planned and ongoing galaxy redshift surveys aim at pre- 
cisely measuring the anisotropic power spectrum and/or 
two-point correlation function in redshift space. Thus, 
the accurate theoretical modeling of anisotropic power 
spectrum is crucial and needs to be developed toward 
future observations. 

The purpose of this paper is to address these issues 
based on the the analytical treatment of non-linear grav- 
itational clustering. In the single-stream limit, cosmo- 
logical evolution of the mass distribution consisting of 
the cold dark matter (CDM) and baryon is described 
by the coupled equations for irrotational and pressure- 
less fluid zz]. Recently, a detailed study on the stan- 
dard treatments of perturbation theory has been made 
[23l425|. a, nd several improved treatments have been pro- 
posed [26|438] , showing that a percent-level accuracy can 
be achieved for the predictions of power spectrum or two- 
point correlation function in real space. With a help of 
these, in this paper, we will develop a model of redshift 
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distortion, and compute the matter power spectrum in 
redshift space, with a particular attention to the BAOs. 
Also, we discuss the impact of model imcertainty of the 
redshift distortion on the acoustic-scale measurement of 
BAOs and the estimation of growth rate parameter. 

This paper is organized as follows: In Sec. we start 
with writing down the relation between real space and 
redshift space, and derive an exact expression for mat- 
ter power spectrum in redshift space. Wc then consider 
the existing theoretical models of redshift distortion, and 
compare those with N-body simulations in Sec. lIIIl show- 
ing that non-negligible discrepancy appears at the scales 
of BAOs. In Sec. IIVI non-linear corrections relevant to 
describe the small discrepancies are derived based on the 
exact expression of rcdshift-space power spectrum. The 
new model of redshift distortion including the corrections 
reproduces the N-body simulations quite well in both the 
monopole and quadrupole components of redshift-space 
power spectrum. In Sec. |Vl the relevance of this new 
model is discussed in details, especially for measurements 
of acoustic scales and growth rate parameters. The po- 
tential impact of the model of redshift distortion on fu- 
ture constraints on modified gravity and dark energy is 
also estimated based on Fisher matrix formalism. Fi- 
nally, our important findings are summarized in Sec. I VII 

Throughout the paper, we assume a flat ACDM model 
and adopt the fiducial cosmological parameters based on 
the five-year WMAP results = 0.279, Qa = 

0.721, llb/f^m = 0.165, h = 0.701, = 0.96, erg = 0.817. 
In order to compare our analytic results with N-body 
simulations, the data are taken from Ref. [sj, in which 30 
independent N-body simulations of the 512^ particles and 
cubic boxes of side length 1/i^^Gpc were carried out with 
initial conditions created by 2LPT code [i^ at Zinit = 31, 
adopting the same cosmological parameters as mentioned 
above. 



II. POWER SPECTRUM IN REDSHIFT SPACE 

Let us first recall that the redshift distortion arises 
from the apparent mismatch of galaxy position between 
real and redshift spaces caused by the contamination of 
the peculiar velocities in the redshift measurement. For 
distant galaxies, the position in real space, r, is mapped 
to the one in redshift space, s, as 



s = r 



aH{z) 



(1) 



where the unit vector 2 indicates the line-of-sight direc- 
tion, and quantity Vz is the line-of-sight component of 
the velocity field, i.e, Vz — v -z. The quantities a and H 
are the scale factor of the Universe and the Hubble pa- 
rameter, respectively. Then, the density field in redshift 
space, i5*^^^(s), is related to the one in real space, 5{r) 
through the relation {1 -I- 5^'^\s)}(fs = {1 -f 5{r)}(fr, 



which leads to 

<5(s)(s) 



ds 
dr 



{l + 5[r))-l. 



(2) 



The Fourier transform of this is given by 



aH{z) J 

(32 

where the quantity /i is the cosine of the angle between z 
and k. Here, we used the fact that the Jacobian \ds/dr\ 
is written as 1 -I- VzVz/{aH). 

From this, the power spectrum of density in redshift 
space becomes [41 [ 



X {6{r) + fWzUzir)} {S{r') + fWzUz{r')}), (4) 

where x = r r' and (• • • ) is the ensemble average. 
We defined Uz{r) = —Vz{r)/{aHf) and Auz = Uz{r) — 
Uz{r'). The function / is the logarithmic derivative of lin- 
ear growth function D{z) given by / = d\xiD[z) / d\na. 
This is the exact expression for power spectrum in red- 
shift space, and no dynamical information for velocity 
and density fields, i.e., Euler equation and/or continuity 
equation, is invoked in deriving this equation. 

In the expression (|4]), the power spectrum is written 
as function of k and and is related to the statistical 
average of real-space quantities in a complicated manner, 
but qualitative effects on clustering amplitude of power 
spectrum are rather clear, i.e., enhancement and damp- 
ing, well known as Kaiser and Finger-of-God effects. The 
Kaiser effect basically comes from the braces in the right 
hand side of the expression (|4]) , which represents the co- 
herent distortion by the peculiar velocity along the line- 
of-sight direction. In linear theory, the relation u = 5 
holds and the strength of clustering anisotropics is con- 
trolled by the growth rate parameter /. This is the basic 
reason why the redshift distortion attracts much atten- 
tion as a powerful indicator for growth of structure. On 
the other hand, Finger-of-God effect roughly comes from 
the factor e-'-ki^f^u^ ^j^g random- 

ness of peculiar velocities, de-phasing arises and it leads 
to the suppression of clustering amplitude. The appar- 
ent reduction of amplitude becomes especially significant 
around the halo forming regions. 

Of course, these two effects cannot be separately 
treated in principle, and a mixture of Kaiser and Finger- 
of-God effects is expected to be significant on trans-linear 
regime, where a tight correlation between velocity and 
density fields still remains. This is of particular impor- 
tance for the accurate modeling of BAOs. Before address- 
ing detailed modeling, however, we will first consider cur- 
rently existing models of redshift distortion, and examine 
how these models fail to reproduce the major trends of 
BAO features in redshift space. 
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III. EXISTING MODELS OF REDSHIFT 
DISTORTION 

A. Perturbation theory description 

Let us first examine the pcrtmbation theory (PT) 
based model of redshift distortion. We here specifically 
deal with the two representative models: one-loop PT 
calculations for redshift-space power spectrum from stan- 
dard PT and Lagrangian PT. 

The standard PT usually implies a straightforward ex- 
pansion of the the cosmic fluid equations around their lin- 
ear solution, assuming that the amplitudes of density and 
velocity fields are small. This treatment is also applied to 
the evaluation of redshift-space power spectrum (|H) , and 
the resultant expressions for one-lo op p ower spectrum is 
schematically summarized as (see j29l . l42j for complete 
expressions) 



d(S) 
SPT 



{k, = (1 + //i^)^ Pii„(fc) + P^^Lpik, /i), (5) 



The first term in the right-hand side is the linear-order 
result of the redshift-space power spectrum, and the fac- 
tor (1 + /^^)^ multiplied by the linear power spectrum 
Piin indicates the enhancement due to the Kaiser effect. 
The second term Pi-i^ap represents a collection of the 
leading-order mode-coupling terms called one-loop cor- 
rection, arising both from the gravitational clustering 
and the redshift distortion. This term is basically of the 
forth order in linear-order density or velocity fields, and is 
roughly proportional to PunA^ with A'^ = A:''Pii,j/(27r^). 

On the other hand, the Lagrangian PT description of 
the redshift-space power spectrum is obtained in some- 
what different way. Intuitively, we rewrite the exact ex- 
pression ^ in terms of the displacement vector, and 
the perturbative expansion is applied to the displace- 
ment vector. Although a naive perturbative treatment 
merely reproduces the standard PT result ([5]), Ref. [l^ 
has applied a partial expansion, and some of the terms 
has been kept in some non-perturbative ways. The re- 
sultant expressions for power spectrum in redshift space 
becomes 



lpt(^' /^) 



X 



^s?t(^> + (1 + /m')' {1 + /(/ + 2)m'} k'al,. 



(6) 

where the quantity cr^ j;„ is the linear-order estimate of 
the one-dimensional velocity dispersion given by 



v,lin 



(7) 



The exponential prefactor in Eq. ^ can be regarded as 
the result of non-perturbative treatment, and in redshift 
space, this term accounts for the non-linear damping of 
the BAOs arising both from the gravitational clustering 
and Finger-of-God effect of redshift distortion. 



Fig. [T] compares the PT based models of redshift dis- 
tortion with N-body simulations of Ref [Sj]. Left and 
right panels respectively show the monopole = 0) and 
quadrupole {£ — 2) moments of power spectrum divided 
by the smooth reference spectrum at different redshifts, 
z = 3, 1 and 0.5 (from top to bottom). The reference 
spectrum P^^no-wiggiG(^) calculated from the no- wiggle 
approximation of the linear transfer function in Ref. psj , 
taking account of the linear-order result of Kaiser effect. 
The multipole moment of two-dimensional power spec- 
trum is defined by 



2^+1 



diiP^^\k,ix)V,{ii), (8) 



with Vi{fJ.) being the Legendre polynomials. 

As it has been repeatedly stated in the literature 
[H m, m. Hi], the standard PT treatment is not suf- 
ficiently accurate to describe the BAOs. Fig. [T] confirms 
that this is indeed true not only in real space, but also 
in redshift space. While the power spectrum amplitude 
of N-body simulations tends to be smaller than that of 
the linear theory prediction (dotted), the predicted am- 
plitude of standard PT generally overestimates the N- 
body results, and it exceeds the linear prediction on small 
scales. Compared to the results in real space, the dis- 
crepancy between prediction and simulation seems a bit 
large. Contrastingly, in the Lagrangian PT calculation, 
the amplitude of power spectrum is rather suppressed, 
and a better agreement between prediction and simula- 
tion is achieved at low-fc. This is due to the exponen- 
tial prefactor in Eq. As a trade-off, however, the 
predicted amplitude at higher k modes largely underesti- 
mates the result of N-body simulations. Further, a closer 
look at first peak of BAOs around k ~ 0.05 — 0.1 ZiMpc^^ 
reveals a small discrepancy, which becomes significant as 
decreasing the redshift and can produce few % errors in 
power spectrum amplitude. 

These results indicate that the existing PT based ap- 
proaches fail to describe the two competitive effects of 
redshift distortion in the power spcctrum[65]. A proper 
account of these is thus essential in accurately modeling 
BAOs. 



B. Phenomenological model description 

Next consider the phenomenological models of redshift 
distortion, which have been originally introduced to ex- 
plain the observed power spectrum on small scales. Al- 
though the relation between the model and exact expres- 
sion (|4]) is less clear, for most of the models frequently 
used in the literature, the redshift-space power spectrum 
is expressed in the form as (e.g., |4ll. |45| - [50| ) 



P'^^^(fc, /i) = DFoG[knf<Jv] PKaiscr(fc,A*): 



(9) 



where the term PKaisci(fc, A*) represents the Kaiser effect, 
and the term -DfoG [k fJ- f ctv] indicates the damping func- 
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FIG. 1: Ratio of power spectra to smoothed reference spectra in redshift space, P^'^' (fc) / Pi^^a-wiggid^) ■ N-body results are taken 

from the wmap5 simulations of Ref. . The reference spectrum P^ no_„iggi(, is calculated from the no- wiggle approximation of 
the linear transfer function, and the linear theory of the Kaiser effect is taken into account. Short dashed and dot-dashed lines 
respectively indicate the results of one-loop PT and Lagrangian PT calculations for redshift-space power spectrum (Eqs. ((5]) 
and ®). 
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FIG. 2: Same as in Fig. [T] but we here plot the results of phenomenological model predictions. The three different predictions 
depicted as solid, dashed, dot-dashed lines are based on the phenomenological model of redshift distortion (|9]) with various 
choices of Kaiser and Finger-of-God terms (Eas. (|10p and (|lip ). Left panel shows the monopole power spectra = 0), and 
the right panel shows the quadrupole spectra {£ = 2). In all cases, one-dimensional velocity dispersion CTv was determined by 
fitting the predictions to the N-body simulations. In each panel, vertical arrow indicates the maximum wavenumber ki% for 
improved PT prediction including up to the second-order Born approximation (see Eq. H12p for definition). 
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FIG. 3: Redshift evolution of velocity dispersion o"v deter- 
mined by fitting the predictions of monopole and quadrupole 
power spectra to the N-body results. While the solid lines rep- 
resent the linear theory prediction, the symbols indicate the 
results obtained by fitting models of redshift distortion with 
various choices of Kaiser and damping terms (see Fig. [2J . 



tion which mimics the Finger-of-God effect. The quan- 
tity cTv is the one-dimensional velocity dispersion defined 
by fly ~ {u1{0)). The variety of the functional form for 

-^Kaiser 

(k,^) and £>FoG[fc M / fv] are summarized as fol- 
lows. 

The Kaiser effect has been first recognized from the 
linear-order calculations [sij . from which the enhance- 
ment factor (1 -I- //i^)^ is obtained (see Eq. ([S])). As a 
simple description for the Kaiser effect, one may naively 
multiply the non-linear matter power spectrum by this 
factor, just by hand. Recently, proper account of the 
non-linear effect has been discussed [4l|, |45|, and non- 
linear model of Kaiser effect has been proposed using the 
real-space power spectra. Thus, we have 



Pk aiscr 



linear 



Pssik) + 2ffj,^Pseik) + p^i^Peeik) ; non-linear 

(10) 

Here, the spectra Pss, Pee, and Pse denote the auto power 
spectra of density and velocity divergence, and their cross 
power spectrum, respectively. The velocity divergence 9 
is defined by 6l = Vm -\/v/{aHf) [66|. 

On the other hand, the functional form of the damp- 
ing term can be basically modeled from the distribution 
function of one-dimensional velocity. Historically, it is 
characterized by Gaussian or exponential function (e.g.. 



[47H50|), which lead to 

exp(— x^) ; Gaussian 



1/(1 + ; Lorentzian 



(11) 



Note that there is analogous expression for exponential 
distribution, i.e., Z^fogN = l/(l+a;^/2)^ but the re- 
sultant power spectrum is quite similar to the one adopt- 
ing the Lorentzian form for the range of our interest, 
X < 1. Since the Finger-of-God effect is thought to be a 
fully non-linear effect, which mostly comes from the viri- 
alized random motion of the mass (or galaxy) residing at 
a halo, the prediction of (Tv seems rather difficult. Our 
primary purpose is to model the shape and structure of 
acoustic feature in the power spectrum, and the precise 
form of the damping is basically irrelevant. We thus re- 
gard CTv as a free parameter, and determine it by fitting 
the predictions to the simulations or observations. 

Fig. [5] compares the phenomenological models of red- 
shift distortion with combination of Eqs. (jTU]) and pT|) 
with N-body simulations. In computing the redshift- 
space power spectrum from the phenomenological mod- 
els, we adopt the improved FT treatment by Refs. (ssl . 
H^l, and the analytic results including the corrections up 
to the second-order Born approximation are used to ob- 
tain the three different power spectra Pss, Pse and Pgg. 
The accuracy of the improved FT treatment has been 
checked in details by Refs. [s^], and it has been shown 
that the predictions of Pss reproduce the N-body results 
quite well within 1% accuracy below the wavenumber 
fci%, indicated by the vertical arrows in Fig. [51 This has 
been calibrated from a proper comparison between N- 
body and FT results and is empirically characterized by 
solving the following equation t25|, iMl : 



(12) 



with C = 0.7 and being linear matter spectrum. 
Note that the 1% accuracy of the improved FT prediction 
at z = 3 has reached at A: ~ 0A7hM-pc~^ , outside the plot 
range. We basically use this criterion to determine CTv, 
and fit the predictions of both monopole and quadrupole 
spectra to the N-body results in the range < fc < ki%. 

Since we allow CTv to vary as a free parameter, the 
overall behaviors of the model predictions reproduce with 
N-body results, and the differences between model pre- 
dictions are basically small compared to the results in 
the FT description. However, there still exist small but 
non-negligible discrepancies between N-body results and 
model predictions, which are statistically significant, and 
are comparable or exceed the expected errors in upcom- 
ing B AO measurements 34] . Although the agreement is 
somehow improved when we adopt the non-linear model 
of ^Kaiser, there still remains discrepancies of few % in 
monopole and 5 % in quadrupole moments of power spec- 
trum amplitudes. These are irrespective of the choice of 
damping function -DfoG- 



6 



Furthermore, the fitted results of CTv show somewhat 
pecuHar behavior. Fig. [3] plots the fitted values of as 
function of redshift (symbols) , which significantly deviate 
from linear theory prediction (solid line) as increasing the 
redshifts. This is in contrast with a naive expectation, 
and indicates that the model based on the expression (|9]) 
misses something important, and needs to be reconsid- 
ered. 



IV. IMPROVED MODEL PREDICTION 
A. Derivation 

Comparison in previous section reveals that even in the 
models with fitting parameter, a small but non-negligible 
discrepancy appears at the scales of BAOs, where the 
choice of the damping function I3fog[2^] does not sensi- 
tively affect the predictions. This implies that there ex- 
ists missing terms arising from the non-linear mode cou- 
pling between density and velocity fields, and those cor- 
rections alter the acoustic feature in redshift-space power 
spectrum. In this section, starting with the exact expres- 
sion (j4]), we derive non-linear corrections, which are rel- 
evant to explain the modulation of acoustic features in 
redshift space. 

First recall that the expression Q is written in the 
form as 

p(s)(A:,^) = f d^a-e^fc-a^^eJi^i^^As), (13) 



where the quantities ji , Ai (i 
given by 



1,2,3) are respectively 



ji = -i knf, 
M = Uz(r) - W2(r'), 
^2 -<5(r) + /V,u,(r), 
^3 -'5(r') + /V.7/.(r'). 

We shall rewrite the ensemble average {e^^^^ A-^A^i) in 
terms of the cumulants. To do this, we use the relation 
between the cumulant and moment generating functions. 
For the stochastic vector field A = {Ai,A2, A^}, we have 
(e.g., MM)- 



) = exp 



(14) 



with j being arbitrary constant vector, j = {ji,j2, js}. 
Taking the derivative twice with respect to j2 and j3, and 
we then set ^2 = = 0. We obtain 41 1 



(e^'i^M2^3) =exp{(e^'^^Oc} 

X [(e^'^^M2^3)c + {e^'^'A2)c{e''^'A^), 
Substituting this into Eq. (|13p . we arrive at 



(15) 



ikx 



exp{(e^^^Oc} 



This expression clearly reveals the coupling between 
density and velocity fields associated with Kaiser and 
Finger-of-God effects. In addition to the prefactor 
exp {(e^^"^Oc}j the ensemble averages over the quantities 
A2 and A3 responsible for the Kaiser effect all includes 
the exponential factor e^^'^^, which can produce a non- 
negligible correlation between density and velocity. 

Comparing Eq. (|16|) with the expression ([9]) with (fTO|) 
and (fTT|) . we deduce that the phenomenological models 
discussed in Sec. lIII 5] miss something important, and are 
derived based on several assumptions or treatments: 



• In the integrand of Eq. (|16p , while taking the limit 
ji ^ in the bracket, we keep ji 7^ in the expo- 



X [{e^'^'A2A3)e + {e''^'A2)c{e''^'A3)c 



(16) 



nent of the prefactor. 



• For cumulants (v4")c = {[uzir) - Uz(r')]")c of 
any integer value n, the spatial correlations be- 
tween different positions are ignored, and the non- 
vanishing cumulants are assumed to be expressed 
as {Ai)c — 2(m")c — 2cncr''^ for even number of n, 
with c„ being constants. 



• To further obtain the Gaussian or Lorentzian forms 
of the damping function DpocN, '^^ assume that 
the conditions, c„ = except for C2 = 1, or, C2n = 
(2n — 1)! and C2n-i = 0, are fulfilled. 



In the above, the last two conditions play a role for 
specifying the damping function, and they mainly affect 
the broadband shape of the power spectrum. On the 
other hand, the first condition leads to the expression of 
^'Kaiscr(^), which Can add the most dominant contribu- 
tion to the acoustic feature in power spectrum. Since 
the choice of the damping function (jlll) is presumably a 
minor source for discrepancies between the model pre- 
dictions and simulations, taking the limit ji — > in 
the bracket would be the main reason for discrepancy. 
In this respect, the terms involving the exponential fac- 
tor can produce additional contributions to the spectrum 
PKaiscr(^), which are responsible for explaining the mod- 
ulated acoustic peak and trough structure in redshift 
space. 

Let us now derive the corrections to PKaisGi (^). To do 
this, we keep the last two conditions, and perturbatively 
treat the terms inside the bracket of Eq. (|16p . This treat- 
ment is reasonable, because the modification of acoustic 
features should be small for relevant scales of BAOs. On 
the other hand, the factor exp{(e^^^i)c} is most likely 
affected by the virialized random motion of the mass 
around halos, and seems difficult to treat it perturba- 
tively. Here, regarding the quantity ji as a small expan- 
sion parameter, we perturbatively expand the terms in 
the bracket of the integrand. Up to the second order in 
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ji, we have 

~ {A2A3,) + ji{AiA2A-i)c 

+ {AIA2A:,), + {A^A2)MiM)c] + 0{]l). 



(17) 



In the above, the term {A\A2A^) c turns out to be higher 
order when we exphcitly compute it employing the per- 
turbation theory calculation, and is roughly proportional 
to ©(PjfjJ. We thus drop the higher-order contribution, 
and collect the leading and next-to-leading order contri- 
butions. Then, Eq. can be recast as 

P(=^) {k,^i) = Df,g [k^i f a,] {P^a ik) + 2f Pse [k) 
+ ffi^Pggik)+Aik,^i) + B{k,n)y (18) 

Here, we replaced the exponential prefactor 
exp{(e-'^"^i)c} with the damping function UfoG- 
The corrections A and B are respectively given by 



Bik,fi) = jf 



{AiA2)c{AiA3)c 



In terms of the basic quantities of density 6 and velocity 
divergence 9 = —\/v/(aHf), they are rewritten as 



A{k,|,) = {k^lf) 



(27r)3 p2 

X {P,(p, k p, -k) - B,{p, k, -k - p)} , (19) 
2 f d?P 



B{k,f,) = (fc^ fr J AP_F{p)F{k - p) ■ 



(20) 



where the function B^ is the cross bispectra defined by 



e(fci) { 5{k2) + / -^0{k2)^ |<S(fc3) + / -^^(fca) 

{2TT)HD{ki + k2 + k;) P,(fci, fc2, fes). (21) 



In deriving the expression (|18p , while we employed the 
low-fc expansion, we do not assume that the terms Ai 
themselves are entirely small. In this sense, the expres- 
sions (HH]) , (HH) and (pn|) still have some non-perturbative 
properties, although the new corrections A and B ne- 
glected in the previous phenomenological models are ex- 
pected to be small, and can be treated perturbatively. 
In Appendix \^ based on the standard PT treatment, 
we summarize the perturbative expressions for the cor- 
rections and (PO)) . in which the three-dimensional 
integrals are reduced to the sum of the one- and two- 
dimensional integrals. 
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FIG. 4: Contributions of power spectrum corrections com- 
ing from the A and B terms divided by the smooth refer- 
ence power spectrum, P^lijk) / P^%_^,^^Jk) (Eq. m)- We 
adopt the Gaussian form of the damping function DfoG with 
(Tv computed from linear theory (see Eq.((7])). Left and right 
panels respectively show the monopole and quadrupole power 
spectra at redshifts z = 3 and 1. 



To see the significance of the newly derived terms A 
and B, we evaluate the monopole and quadrupole con- 
tributions to the functions defined by 



(22) 



The results are then plotted in Fig. |31 divided by the 

(3) 

smoothed reference spectrum, no-wiggie(^)- plotting 
the results, we specifically assume the Gaussian form of 
DfoG: s-nd adopt the linear theory to estimate (Tv (see 
Eq. dZl)). 

The corrections coming from the A term show oscilla- 
tory behaviors, and tend to have a larger amplitude than 
those from the B term. While the corrections from the 
B term are basically smooth and small, they still yield 
a non-negligible contribution, especially for quadrupole 
power spectrum. Although the actual contributions of 
these corrections to the total power spectrum are deter- 
mined by the fitting parameter CTv, and thus the resul- 
tant amplitudes shown in Fig. 2] do not simply refiect 
the correct amplitudes, the new corrections A and B can 
definitely give an important contribution to the acoustic 
feature in power spectrum. 

Finally, it is interesting to note that while the new for- 
mula for redshift-space power spectrum IjlSp would be 
applicable to the non-linear regime where the standard 
PT calculation breaks down, the resultant expression it- 
self is similar to the one for redshift-space power spec- 
trum in the one-loop standard PT. The one-loop power 
spectrum in redshift space, Pgprj,(A:, /i) given at Eq. (O, 
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FIG. 5: Same as in Fig. [21 but we here adopt new model of redshift distortion ((TSJ. Solid and dashed lines represent the 
predictions for which the spectra P^a, P^b and Pee are obtained from the improved PT including the correction up to the 
second-order Born correction, and one-loop calculations of the standard PT, respectively. In both cases, the corrections A and 
B given in Eqs. (I19|l and (|20|l are calculated from standard PT results (see Appendix The vertical arrows indicate the 
maximum wavenumber k\% defined in Eq. ([12}, for standard PT and improved PT (from left to right). 



can be formally recast as 

P^%{k,fi) = {l- (fc^Mjin)'} {Pss{k)+2ffi^Pse{k) 
+ffi^P9e{k)} + A{k, ix) + B{k, fi) + C{k, fi). (23) 

Note that each term in the above expression should be 
consistently evaluated using the perturbative solutions 
up to the third order in d and 0, and as a result, only 
the leading-order corrections just proportional to PunA^ 
(or equivalently the forth order in J*-^-*) are included in 
the one-loop power spectrum. Here, the function C is 
defined by 

C{k,f.) = {kf,fr I f0^S^{k-p-q)^Pee{p) 
X { Pss {q)+2f ^2 p,g (q) + f ^4 p^^ I 

X Pii„(p)Pii„(<z) (24) 

with fip ^ Pz/\p\ and /i^ = qz/\q\. The second equality is 
valid for the one-loop PT calculation. Hence, if we adopt 
either of Lorentzian or Gaussian form in Eq. (jlip and just 
expand it in powers of its argument, the new formula (fT8|) 
reduces to the one-loop result (1^5)) just dropping the term 
C. 

The C term is originated from the spatial correlation 
of the velocity field, and is obtained through the low-A: 



expansion of the exponential prefactor exp{(e^^"*^)c} in 
Eq. ([T5)) . For the scales of BAOs, the C term monoton- 
ically increases the amplitude of power spectrum, and 
it does not alter the acoustic structure drastically. In- 
deed, our several examinations reveal that the effect of 
this can be effectively absorbed into the damping func- 
tion D[fc/i/(Tv] with varying the velocity dispersion tJv. 
Rather, the main drawback of the standard PT expres- 
sion (|23p comes from a naive expansion of all the terms 
in the exact formula (|3]), which fails to describe the del- 
icate balance between the Finger-of-God damping and 
the enhancement from Kaiser effect and non-linear grav- 
itational growth. As we will see in next subsection, both 
keeping the damping term DpoG and including the cor- 
rections A and B seem essential, and with this treatment, 
even the standard PT calculation of the power spectrum 
can give a excellent result which reproduces the N-body 
simulations fairly well. 



B. Comparison with N-body simulations 

We now compare the new prediction of redshift-space 
power spectra with the result of N-body simulations. 
Fig. [5] shows the monopole (left) and quadrupole (right) 
power spectra divided by their smooth reference spectra. 
The analytical predictions based on the model are 
plotted adopting the Gaussian form of the Finger-of-God 
term DpoQ[kf ^a^], and the velocity dispersion CTv is de- 
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FIG. 6: Contribution of each term in redshift-space power 
spectrum. For monopole {I — 0, left) and quadrupole 
{£ = 2, right) spectra of the improved model prediction at 
z = 1 shown in solid lines of Fig. (5] we divide the to- 
tal power spectrum P^^^^^ (solid) into the three pieces as 



p(S) _ p' 
total Kaiser 



3(S) 



I p(S) , p(S) 
' corr,A corr,. 



and each contribution is 



separately plotted dividing by smoothed reference spectra, 



3(S) 
!,no-wigglG ■ 



Here, the spectrum ^Kaiser (dotted) is the con- 
tribution of non-linear Kaiser term pO|) convolved with the 



Finger-of-God damping DfoG, and the corrections 



(S) 



and 



3(S) 



are those given by Eq. (p2l 



termined by fitting the predictions to the N-body results. 
In computing the predictions, the A and B terms are 
calculated from the one-loop standard PT results in Ap- 
pendix]^ while the spectra Pss, Pse and Pgg are obtained 
from improved PT in solid lines, and from standard PT 
in dashed lines. 

Compared to Figs. [Hand [21 the agreement between N- 
body simulations and predictions depicted as solid lines 
becomes clearly improved, and the prediction includ- 
ing the corrections faithfully traces the N-body trends 
of acoustic feature, especially around k = 0.05 ~ 
0.15/iMpc~^, where the phenomenological model shows 
a few % level discrepancy. A remarkable point is that a 
reasonable agreement basically holds over the range be- 
low the critical wavenumber ki% calibrated in real space 
(vertical arrows, Eq. (|12p for definition). This is also 
true for the case adopting one-loop standard PT to com- 
pute Pss^ Pse and Pqq (dashed lines), and the range of 
agreement is wider than that of the existing PT-based 
models in Sec. IIII Al 

In Fig. [6l to see the significance of the contributions 
from corrections A and B, we divide the improved PT 
prediction of power spectra P'^^\k) at z = 1 into the 



5(S) 



5(S) 



and P 



(S) 



which are 



three pieces as P^^;,,,,, . ,„„_a - corr.B' 
separately plotted as dotted, long-dashed, and short 

(s) 

dashed lines, respectively. The power spectrum PKaisor 
is the contribution of the non-linear Kaiser term given 
in Eq. (|10p. convolved with the damping function UfoG- 
The spectra P^orr A ^'^^ -^corr B represent the actual con- 
tributions of the corrections A and B defined by Eq. 



standard PT 
Improved PT 




Linear 



FIG. 7: Same as in Fig. [S] but we here adopt new model 
of redshift distortion in estimating ctv. The filled triangle 
and circles are the results obtained from predictions based on 
standard PT and improved PT calculations, respectively (see 
dashed and solid lines in Fig. [5] 



with fitted value of CTv The corrections A and B give 
different contributions in the amplitude of monopole 
and quadrupole spectra, and their total contribution can 
reach ^ 10% and ^ 40% for monopole and quadrupole 
spectra at fc < 0.2/iMpc~^, respectively. Thus, even 
though the resultant shape of the total spectrum P'^^^ (k) 
apparently resembles the one obtained from phenomeno- 
logical model, the actual contribution of the corrections 
A and B would be large and cannot be neglected. 

Note, however, that a closer look at low-z behavior re- 
veals a slight discrepancy around fc ^ 0.15/iMpc~^ and 
0.22/iMpc~^ in the monopole spectrum. Also, discrep- 
ancies in the quadrupole spectrum seems bit large, and 
eventually reach ~ 5% error in some wavenumbers at 
z = 0.5. This is partially ascribed to our heterogeneous 
treatment on the corrections A and B using the standard 
PT calculations. It is known that the standard PT result 
generically gives rise to a strong damping in the BAOs, 
and it incorrectly leads to a phase reversal of the BAOs. 
Thus, beyond the validity regime of the standard PT, the 
predictions including the small corrections tend to over- 
smear the acoustic feature, leading to a small discrepancy 
shown in Fig. [S] 

Another source for the discrepancies may come from 
the effect of finite-mode sampling caused by the finite 
boxsize of the N-body simulations. As advocated by 
Refs. [25l. [52I ■ due to the finite number of Fourier modes, 
the matter power spectrum measured from N-body sim- 
ulations may not agree well with the predictions of linear 
theory nor standard PT even at very large scales, and 
tends to systematically deviate from them. While we 
follow and extend the procedure of Ref.fl^ to correct 
this effect in redshift space, it relies on the leading-order 
calculations of standard PT, and the correction for finite- 
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mode sampling has been restricted to the low-fc modes, 
k < 0.1/iMpc~i [3J]. Hence, the high-fc modes of the 
power spectrum plotted here may be affected by the ef- 
fect of finite-mode sampling, and it would be significant 
for higher-multipole spectrum because of its small am- 
plitude. This might be still serious even with the 30 
independent data of N-body simulations. 

Perhaps, the best way to remedy these discrepancies 
at low-z is both to apply the improved PT treatment 
to the corrections A and B, and to consider the higher- 
order contributions for correcting the effect of finite-mode 
sampling over the relevant range of BAOs. The complete 
analysis along the line of this need some progress and is 
beyond the scope of this paper. Nevertheless, it should 
be stressed that the model given by Eq. captures 
several important aspects of redshift distortion, and even 
the present treatment with standard PT calculations of 
the corrections A and B can provide a better description 
for power spectra. In Fig. [3 we plot the fitted values 
of the velocity dispersion obtained from the new predic- 
tions shown in Fig. [Sj The redshift dependence of the 
fitted results roughly matches physical intuition, and is 
rather consistent with the linear theory prediction. This 
is contrasted to the cases neglecting the corrections (see 
Fig.©. 

As another significance, we plot in Fig. [5] the 
quadrupole-to-monopole ratios for redshift-space power 
spectra. The new model predictions using standard and 
improved PT calculations (solid and dashed) are com- 
pared with those neglecting the corrections A and B (dot- 

dashed). The amplitude of the ratio V^o basically 
refiects the strength of the clustering anisotropics, and is 
proportional to (4//3 + 4/77)/(l + 2//3 + p /b) in the 
limit fc — (e.g., [llH, [5l[). One noticeable point is that 
the N-body results of quadrupole-to-monopole ratio do 
exhibit an oscillatory behavior, and the model including 
the corrections reproduces the N-body trends fairly 
well. On the other hand, the phenomenological model 
neglecting the corrections generally predicts the smooth 
scale-dependence of the ratio / and thus it fails 
to reproduce the oscillatory feature. Since this oscillation 
is originated from the acoustic feature in BAOs, Fig. [8] 
implies that the quadrupole-to-monopole ratio possesses 
helpful information not only to constrain the growth-rate 
parameter /, but also to determine the acoustic scales. In 
other words, any theoretical template for redshift-space 
power spectrum neglecting the corrections A and B may 
produce a systematic bias in determining the growth-rate 
parameter /(z), Hubble parameter H{z) and angular di- 
ameter distance Da{z), which we will discuss in details 
in next section. 



V. IMPLICATIONS 

The primary science goal of future galaxy surveys is to 
clarify the nature of late-time cosmic acceleration, and 
thereby constraining the parameters Da{z), H{z) and 




k [h Mpc-i] 

FIG. 8: Quadrupole-to-monopole ratios for redshift-space 
power spectrum, P2^\k)/ PQ^\k), given at z — 3, 2, 1, and 
0.5 (from top to bottom). Solid and dashed lines respectively 
represent the predictions based on new model of redshift dis- 
tortion combining improved PT and standard PT calculation 
to estimate the three different power spectra Pss , Pse and Pee . 
Dot-dashed lines are the results based on the phenomeno- 
logical model neglecting the corrections, which correspond to 
solid lines in Fig.[2](i.e., non-linear Pxaiser -|- Gaussian Dfog)- 
The vertical arrows indicate the maximum wavenumber ki% 
for standard PT (left) and improved PT (right). 

f{z) through a precise measurement of BAOs in redshift 
space would be the most important task. However, these 
constraints may be biased if we use the incorrect model 
of redshift distortion as theoretical template fitting to 
observations. In this section, we explore the potential 
impact on the uncertainty and bias in the parameter es- 
timation for Da{z), H{z) and f{z). 



A. Recovery of parameters Da, H and / 

Let us first examine the parameter estimation using 
the new model of redshift distortion. Fitting the theo- 
retical template of power spectrum to the N-body data, 
we will check if the best-fit parameters for Da{z), H{z) 
and f{z) can be correctly recovered from the monopole 
and quadrupole moments of anisotropic BAOs. 

We model the power spectrum of N-body simulations 
by 

where the comoving wavenumber k and the directional 
cosine fi for the underlying cosmological model are re- 
lated to the true ones q and v by the Alcock-Paczynski 
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effect through (e.g., HSglSC 



q = k 




^1 



1/2 



(26) 



H 




D 



1 -1/2 



A, fid \ I 2 
1^' 



(27) 



The quantities I?yi,fid and -fffid are the fiducial values 
of the angular diameter distance and Hubble parameter 
adopted in the N-body simulations. For a given set of 
cosmological parameters, the redshift-space power spec- 
trum P^^^f is calculated from Eq. p^ . but wc here treat 
the quantity / as free parameter in addition to the ve- 
locity dispersion a^- Further, to mimic a practical data 
analysis using galaxy power spectrum, we introduce the 
bias parameter &, assuming the linear deterministic rela- 
tion, i.e., dsim — bSm ,67]. Then, fitting the monopole 
and quadrupole power spectra of Eq. to those of the 
N-body simulation at z = 1, we determine the best- fit 
values of Da, H and /, just marginalized over the pa- 
rameters (Tv and b. To do this, we use the Markov chain 
Monte Carlo (MCMC) technique described by Ref. [H], 
and adopt the likelihood function given by 

-21n£ = ^ E {^f2n(fc«)-^fidci(fcn)} 

n £,£'=0,2 

X Covllik^) {p/^L(fcn) - ^Onodcl(fcn)} . 

(28) 

where the quantity Covij' represents the covariance ma- 
trix between different multipoles. The range of wavenum- 
ber used in the likelihood analysis was chosen as fc < 
kmax — 0.205/iMpc~^, so as to satisfy fcmax < As 
for the covariance, we simply ignore the non-Gaussian 
contribution (see Ref. [13] for validity of this treatment), 
and use the linear theory to estimate the diagonal com- 
ponents of the covariance, Cov^^^' , including the effect of 
shot-noise contribution assuming the galaxy number den- 
sity rTg = 5 X 10~^/i^Mpc~^. The explicit expression for 
the covariance is presented in Appendix O We checked 
that the linear theory estimate reasonably reproduces the 
N-body results of the covariance matrix for the range of 
our interest k < 0.3/iMpc^^ at z = I. 

Fig. |9] summarizes the result of the MCMC analy- 
sis assuming an idealistically large survey with Vs = 
20h^^Gpc^ . The two-dimensional contour of the 1-a 
marginalized errors are shown for D^/D^i^d vs H/Hfn 
(bottom left), DA/DAfid vs / (middle left), and / vs 
Da/ Da fid (bottom center). Also, the marginalized pos- 
terior distribution for each parameter are plotted in the 
top left, middle center, and bottom right panels. In each 
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FIG. 9: Results of MCMC analysis using the model of red- 
shift distortion with and without corrections (depicted as blue 
and red lines, respectively). Based on the power spectrum 
template (|25p . we derive the posterior distribution for the pa- 
rameters Da, H and / from the monopole and quadrupole 
spectra of N-body simulations at 2 = 1, marginalized over 
the one-dimensional velocity dispersion CTv and linear bias pa- 
rameter h. Top left, middle center and bottom right show 
the marginalized posterior distribution for Da/ DA,Sid, H/Hud 
and /. Shaded regions indicate the 1% interval around 
the fiducial values. Middle left, bottom left, and bottom 
center plot the two-dimensional 1-a errors on the surfaces 
(H/HfiiJ), {DA/DA,Rd,H/HRi), and (/,i////fid). Note that 
in estimating likelihood function (|28|) . we adopted the lin- 
ear theory to calculate the covariance matrix Cov^^^/ , includ- 
ing the shot-noise contribution with Ug — 5 x 10~*/i"^Mpc~'^ 
and assuming an idealistically large survey volume Vs = 
20/i~^Mpc^ (see Appendix [Cl for explicit expression). 



panel, blue and red lines respectively represent the results 
using the model of redshift distortion with and without 
the terms A and B. 

As it is clear from Fig. [SI the model including the cor- 
rections shows a better performance. Within the 1-a 
errors, which roughly correspond to the precision of a 
percent-level, it correctly reproduces the fiducial values 
of the parameters (indicated by crosses). On the other 
hand, the two-dimensional errors of the results neglecting 
the corrections show a clear evidence for the systematic 
bias on the best-fit parameters. Accordingly, the resul- 
tant value of around the best-fit parameters, given by 
= — 2 ln£, is larger than that of the case including the 
corrections: = 10-1 and 22.2 for the cases with and 
without corrections, respectively. Although the deviation 
from the fiducial values seems somewhat small except for 
the growth-rate parameter /, this is solely due to the fact 
that we only use the monopole and quadrupole power 
spectra. It would be generally significant in the analy- 
sis using the full shape of redshift-space power spectrum. 
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for which the statistical errors are greatly reduced, and 
thereby the systematic biases would be prominent. 



B. Impact of redshift distortion on future 
measurements of Da, H and / 

Given the fact that the robust measurement of Da, H 
and / can be made with the new model of redshift dis- 
tortion, we then move to the discussion on the potential 
impact on the future measurements using the full shape 
of the redshift-space power spectrum. Here, for illustra- 
tive purpose, we consider the two surveys around z = 1, 
with volume K = 4 and 20/i~^Gpc^, and quantitatively 
estimate how the wrong model of redshift distortion leads 
to the incorrect measurements of D^, H and /. 

The fundamental basis to estimate the uncertainties 
and systematic biases on model parameters is the Fisher 
matrix formalism. The Fisher matrix for galaxy survey 
is given by 



" (27r)2 



dk 



dp. 



91nPi^)(fc,A.)ainPi^)(fc,^) 



obs V 



obs ' 

dpj 



,(S) 
'g obs 



(29) 



with rlgai being the number density of galaxies, for which 
we specifically set rlgai = 5 x 10~'^ h^Mpc~^ . The min- 
imum wavenumber available for a given survey, fcmin, is 
set to 2Tr/Vs^^ . Here, the observed power spectrum P^^l 



is given by Eq. (j25p . and we allow to include the influence 
of galaxy biasing adopting the deterministic linear rela- 
tion, 5gai = bSm- Then, we have five parameters in total, 
given by pi = {b, /, a^, Da/Dam^ H/Hud}. Fiducial 
values of these parameters are set as b — 2, f — 0.858, 
DA/DA.Rd = 1 and H/Hf^^ — 1. As for the velocity 
dispersion a^, we use the fitted result to N-body simu- 
lations adopting the new model of redshift distortion in 
Sec. lIVBl and set = 395 km s^^. 

Based on the Fisher matrix (P^ . the systematic bias 
for parameter pi caused by incorrectly modeling theo- 
retical power spectrum is estimated from the following 
formula: 



(30) 



where F'~^ is the inverse Fisher matrix evaluated at the 
fiducial parameter set, but is obtained from an incorrect 
model of redshift distortion as a theoretical template of 
redshift-space power spectrum. The vector Sj is given by 



{2n) 



dkP l\, ^-(fe.M) amp— Mfc,M) 

_i P™"g(fc,^) dpj 



ngP™"s(fc,^) 
1 + ngP™"s(fc,^) 



(31) 



The function P™"s(fc^^) jg the theoretical template 
adopting the incorrect model of redshift distortion. The 
systematic differences in the power spectrum amplitude 
are quantified as P^y^{k, n) = P™"8(fc, p) - P*™°(fc, /i), 

where P^'^^° {k , ij,) is the correct template for redshift- 

(s) 

space power spectrum P^i^g, for which we assume the 
new model of redshift distortion including the terms A 
and B (Eq. p^ 'l. Below, we will quantify the magnitude 
of systematic biases if we incorrectly apply the model of 
redshift distortion neglecting the corrections A and B for 
the power spectrum template. 

Fig. [TU] plots the results of the Fisher matrix calcula- 
tions marginalized over the nuisance parameters b and 
CTv The uncertainties and biases for the best- fit val- 
ues of /, Da and H are estimated assuming fcmax = 
0.12/iMpc-i (left) and 0.2/iMpc^i (right), and the re- 
sults are shown for the surveys with Vg = A h~^Gpc^ 
(open) and 20 h'^^Gpc^ (shade). In each panel, two- 
dimensional contours around the crosses and filled tri- 
angles show the expected 1-cr (68% C.L.) errors around 
the best-fit values adopting the model of redshift distor- 
tion with and without the corrections, respectively. The 
differences between best-fit values (crosses and filled tri- 
angles) represent the systematic biases estimated from 
Eq. (j30p . which remain unchanged irrespective of the 
survey volume Vs. Since the size of marginalized uncer- 

— 1/2 

tainties is proportional to V^s , the systematic bias in 
the best-fit parameters become relatively prominent and 
is considered to be a serious problem if we increase the 
survey volume. Note that similar to the result in Fig. |9l 
there exists a tight correlation of the parameters between 
the growth rate parameter / and quantities Da and H. 
This is consistent with the finding by Ref . [11] , indicating 
that the distinguishing dark energy from modified gravity 
needs another observational constraint. 

Fig. [To] implies that phenomenological model of red- 
shift distortion neglecting the corrections can produce 
a large systematic error. As increasing the maximum 
wavenumber /cmax, the bias on the measurements of an- 
gular diameter distance and Hubble parameter reaches 
~ 1 — 2% error, while the best-fit value for the growth 
rate parameter would be seriously biased with ~ 5% 
error. If we conservatively choose a smaller value of 
kmax ^ 0.12/iMpc~^, these systematics could be still 
within the size of statistical error for surveys with typ- 
ical volume of Vg ^ 4ft,~^Gpc^. However, for a survey 
of larger volume with Vs > 20ft,~'^Gpc^, the systematic 
error on the growth rate parameter becomes outside the 
the marginalized uncertainty. If we aggressively choose 
^max ~ 0.2/iMpc~^ in order to reduce statistical uncer- 
tainties, the systematic biases become definitely serious 
issues in all of the parameters /, Da and H for both 
surveys of volume T4 = 4 and 20/i~'^Gpc'^. Hence, cor- 
rectly modeling redshift distortion would be very crucial 
for both stage-Ill and -IV class surveys defined by the 
Dark Energy Task Force [HI] . 
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DA/DA,f,d f oyo,,,, f 

FIG. 10; Expected two-dimensional contours on marginalized errors around the best-fit values of Da/ DA.&d vs H /Had (bottom 
left), / vs H /Hfid (bottom right) and DA/DA,fid vs / (top left) ai z — 1, obtained from the full shape of redshift-space power 
spectrum. The maximum wavenumber for parameter estimation is chosen as fcmax = 0.12 (left) and 0.2ftMpc~^ (right), so as 
to satisfy the condition fcmax < ki% for standard PT and improved PT, respectively. In each panel, open and shaded contours 
indicate the two dimensional errors for the surveys with Vs — i and 20^~"^Gpc^. 



VI. DISCUSSION AND CONCLUSION 

In this paper, we have investigated the power spec- 
trum in redshift space, and presented a new model of red- 
shift distortion, which is particularly suited for modeling 
anisotropic BAOs around k — ^ 0.3 hMpc~^. Contrary 
to the previous phenomenological modes in which the ef- 
fects of Kaiser and Finger-of-God are separately treated 
in a multiplicative way, the new model includes the cor- 
rections coming from the non-linear coupling between ve- 
locity and density fields, which give rise to a slight up- 
lift in the amplitude of monopole and quadrupole power 
spectra. The model predictions can give a good agree- 
ment with results of N-body simulations, and a percent 
level precision is almost achieved. 

Based on the new model of redshift distortion, we pro- 
ceeded to the parameter estimation analysis, and checked 
if the theoretical prediction correctly recovers the cosmo- 
logical information from the monopole and quadrupole 
spectra of N-body simulations. MCMC analysis revealed 
that while the new model of redshift distortion combin- 
ing the improved PT calculation faithfully reproduces the 
fiducial parameters Da, H and / and the precision can 
reach at a percent level, the model neglecting the correc- 
tions {A and B terms) exhibits a slight offset of the best- 
fit values. In order to estimate the potential impact on 
the future measurement, we have further made the Fisher 
matrix analysis using the full shape of power spectrum 
P'^^^(fc, fi), and found that the existing phenomenological 
models of redshift distortion neglecting the corrections 



produce a systematic error on measurements of the angu- 
lar diameter distance and Hubble parameter by 1 ~ 2%, 
and the growth rate parameter by ~ 5%. This would 
become non-negligible for stage-Ill and -IV class surveys 
defined by the Dark Energy Task Force. Correctly mod- 
eling redshift distortion is thus crucial, and the new pre- 
scription of redshift-space power spectrum presented here 
plays an essential role in constraining the dark energy 
and/or modified gravity from anisotropic BAOs. 

Finally, we note several remaining tasks in practical ap- 
plication to the precision measurement of BAOs. One is 
the improved treatment for calculation of the corrections, 
A and B terms, which needs to evaluate the bispectrum 
of density and velocity fields. In doing this, a system- 
atic treatment using multi-point propagator developed 
by Ref. [s^l would be useful and indispensable. Also, the 
effects of the new contributions to the redshift-space clus- 
tering in the presence of the primordial non-Gaussianity 
and the dark sector interaction would be presumably im- 
portant (e.g., [13, m, H^), and should deserve further 
investigation. Of course, the biggest issue is the galaxy 
biasing. Recent numerical and analytical studies claim 
that the scale-dependent and stochastic properties of the 
galaxy bias can change the redshift-space power spec- 
trum, and the potential impact on the determination of 
the growth-rate parameter would be significant [13, HH . 
A realistic modeling of galaxy biasing relevant for the 
scale of BAOs is thus essential, and a further improve- 
ment of the power spectrum template needs to be devel- 
oped. 
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Appendix A: Perturbation theory calculations for 
correction terms 

In this Appendix, we present the perturbative expres- 
sions for the corrections A and B defined in Eqs. 
and (j20[) , which are originated from the coupling between 
Kaiser and Finger-of-God effects. 

Let us first consider the correction A, which involves 
the bispectrum of density and velocity divergence (see 
Eq. (f2T|) ). Using the perturbative solutions up to the 
second-order, the leading-order result of the bispectrum 
becomes 



S<,(fci,fc2,fc3) = (-2/) 



fr2 



pf) [^ + ^f] G2(fe2,fc3)fli„(fc2)Plin(P3) 



4/G2(fci,fc3) } PUh)PUP3) 



1 



|F2(fei,fc3) 

^f) [p^ikiM) + TTf G2{kiM) } PUki)PUP2) 



(Al) 



with F2 and G2 being the second-order perturbation ker- 
nels given by (e.g., p3. [26l [6^) 

^^(^-^^)-7 + ^U + A:rj+H^J ' 

n , A 3 fci • fc2 / A:i fc2 \ 4 / fci • fc2 V 

Note that the bispectrum (jAip possesses the follow- 
ing symmetries: i?cr(fci, ^2, fcs) = -80^(^1,^3,^2) = 
_Bcr(— fci, — A;2, — fca). Then, substituting the expression 
(IA1[) into the definition (IT^ . the correction A can be 



recast schematically in the form as 

Aik,^,) = -kiJiY.j />: Q„.„(fc,p), (A2) 

where the function Qmn is the scalar function of k and 
p. To further perform the angular integral, we use the 
formulae presented in Appendix |B] (Eq. (jB3p ) , which can 
be obtained by utilizing the rotational covariance of the 
integral. After straightforward but lengthy calculation, 
the correction A{k, ^) is finally reduced to the following 
form: 



(2^ 



dr 1 da; I A„j„(r, x) Pii„(fc; z) + Amn{r, x) P\in{kr; z) | 



Piin (fcVl -h r2 - 2rx; z) 
(1 + r2 - 2rxY 

I 



P\i'a{k;z) I dr Qjnnir) Piinikr; z) 
Jo 



(A3) 



where we introduce the quantities r — k/p and x — {k ■ 
p)/k/p. Note again that fj, is the cosine of the angle 
between line-of-sight direction z and the vector fc, i.e., 
^ — (k ■ z)/k. The non-vanishing components of Amm 



Amn, and Omn sre 



+r{-3 + x^ - 12x^)} , 
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A22 =— {r'^x{li - 41a;2) - A{x + Qx^) 
+r{b + %x^ + A2x'^)] , 



0,12 



112r3 



-3(r2-l)4log 



r + 1 



r - 1 



^23 — ^12, 



A33 = —(1 - 7rx + 6a;2) {-2x + r(-l + Zx'^)] , 



for j4„ 



All = y(a; + r - 2rx^)(3r + 7x- IQrx^), 



a22 =33g^ [2^(9 - 185r2 + 159r4 - GSr^) 



+9(r2 - l)3(7r2 + l)log 



r + 1 



r - 1 



= — ~ + 7a; - 10ra;2), 



023 — ^12) 



^22 = {^S'-r^ + - 81x2) + ^2(^ _ 27a;2 + 54a;'')} , 



^23 = — (l-a;')(r-7a; + 6ra;2), 



A33 = — (r - 7x + 6ra;2)(-2a; - r + Srx^), 



for A^„, and 



an = r2r(19 - 24r2 + Qr"* 

84r 

r + 1 



9(r^ - l)-'log 



r- 1 



033 [2^9 - 109r2 + 63r4 - 27r6) 



+9(r2 - Ifiir^ + l)log 



r + 1 



r - 1 



for 



Next consider the correctfons B. This term is already 
of the order C'({Piin(A;)}2), and the non- vanishing contri- 
bution can be estimated without employing the pertur- 
bative calculations. Just applying the formulae (|B3[) in 
Appendix m to Eq. ((20|) . we obtain 



(A4) 



r 



where Pi2{k) = Pse{k) and P22{k) = PeeCfc). The non- 
vanishing coefficients B^^^ are 



B2\ = ^(^^"ir, 



i?ii = -(x2-l), 



B2^2 = ^(-^-l)^ 



Si2 = ^(x2-lf, 



Bli - -(r + 2a;-3rx2), 
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where we write p = kr and k ■ p = k?rx. Thus, as a 
special case with i = j = z, we get the following formula: 



3r^ 



J (C^P'/(^'P) = ^+(M'Q (B2) 



3^2 

^ __(a;2 - 1)2(6 _ 30rx - + idr^x"^), 



BI2 = o {4a;(3 - hx^) + r(3 - SOx^ + SSx"*)} 



= ^ [-8a; + r {-12 + 36a;2 + 12ra;(3 - 5a 



+r2(3-30a;2 + 35x'*)}] , 



with kz — k ^. 

The above procedure can be generalized to the cases 
of integrals involving arbitrary numbers of multiplicative 
vectors. As a useful formula, we here explicitly write 
down the result summing up the integrals of arbitrary 
scalar functions /„ times the power up to the sixth 
order: 



m.n—0 



Bi 



3r 



— (x^ - 1) r-8a; + r{-l2 + eOa;^ 
16 



+20ra;(3 - 7a;2) + 5r2(l - lAx^ + 21a;'')}] , 



^22 = [8a;(-3 + 5a;2) - 6r(3 - 30a;2 + 35a 

+ 6r2a;(15 - IQx^ + 63a;'') 

{5 - 21a;2(5 - 15^2 + 11a;'')}] . 



The expression (jA4p is still non-perturbative in the 
sense that we do not perturbatively treat the power spec- 
tra Pse and Pee in the integrand. For the leading-order 
calculation, we simply apply the linear-theory calculation 
to these quantities, and replace both Pgg and Pee with 
the linear spectrum Pwa- 



Appendix B: Some useful formulae for integrals 



x/ drr2 / dxikr^G^M fm{Kr,x), (B3) 



The non-vanishing coefficients Gnm as functions of k, r 
and X are summarized as follows: 

Goo = 1, 



G02 = -^{x' - 1), 



Go4 - x(a;2 - l)^ 



Go6 = -Y^(a;^-lf 



In this Appendix, we give the integral formulae used 
in Appendix lAl to derive the perturbative expressions for 
the correction A and B. 

Let us first consider the integral of an arbitrary scalar 
function f{k,p) times some vectors over p. A simple ex- 
ample of the integrand are piPj f(k^p), where subscript 
i, j selects a;-, y- or z-direction. The rotationally invari- 
ant properties of the integral implies that the resultant 
form of the integral is given by 



_cPp_ 

(27r)3 



PiPj f{k, p) ^ P Sij + Q hkj 



(Bl) 



irrespective of the functional form of f{k,p). The coeffi- 
cients P and Q are obtained by contracting the above in- 
tegral with Sij and k^kj, and are the functions of fc = |fc|. 
We have 



P 



Q = 



(2^)2 

(2vr)2 



dr r 



dr 



1 ^2 

dxx — (1 
-1 2 ^ 



x^) f{k,r,x), 



dxx — (3a;2 — l) /(A;, r, x), 
1 2 



Gil = X, 



Gi3 = --a;(a;2 - 1), 



G,, = ^x{x'-lf, 



G22 = -(3x2 _ 



G24 = -^(5a;4-6.a;2 + l), 



G26 = ^(7x2-l)(a;2-lf, 
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G33 - -a;(5a;2-3), 
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G35 = --x{7x^ -lOx^ + 3), 



G44 = -(35x^-30x2 + 3), 



For an analytic estimate of the covariance, we assume 
that the power spectrum is simply written as P'-^-' (fc, /i) = 
(1 + (3 ii'^)'^b'^ Piin{k), where b is the linear bias parame- 
ter and (3 is defined by /3 = f/b. Substituting this into 
Eq. (jCll) . we obtain 



= ^(-21x^ + 35.T^ - 15x2 + 1), 



G55 = -(63x'* - 70x2 ^ ^5^1^ 



Gee = TTT (231x6 
16 



315x'' + 105x2 



5). 



Covo,o(fc) 



Nk 



for (^,f)-(0,0), 



(C2) 



Appendix C: Covariance between multipole power 
spectra 

Here, we give the explicit expression for covariance be- 
tween multipole power spectra used in the MCMC anal- 
ysis in Sec. IV Al 

Neglecting the non-Gaussian contribution, the non- 
vanishing part of the covariance only appears at the di- 
agonal components (i.e., correlation between the same 
Fourier modes), which are expressed as (e.g., Ref. [gl. 

Hi) 



COV£_£'(fc) 



2 (2£-hl)(2£'-hl) 
2 



Covo,2(fc) 



^/3+f/3^ + g/3^ + ^/?^j{6^Pn„(fc)} 

^(tf3+yAb'Plin{k) 



rig V 3 7 



(C3) 



for {£,£') = (0,2) or (2,0), and 



C0V2^2(fc) 



220 „ 90 



5 + ^/3 + ::^/32 + 1122^33. 
2r 7 23r 



2075 
1287 



1 r 1^2 r^2 n /,m2 2 / 220 ^ 30„2\ 

^d,Vei,)VA,){p^''\k,,) + l-] , (CI) x{62p,„(fc)} +=^(5 + — /3-f-/32j 



where Nk is the number of Fourier modes within a given 
bin at k, given by Nk = -iTT k'^Ak/{2T:/Vs^^f with Afc 
and Vs being the bin width and survey volume, respec- 
tively. 



x62pii^(fc) + 



(C4) 



for (^,f) = (2,2). 
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